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H , Abstract 

We present the values of the complete set of planar four-point Master Integrals 
needed for massive Bhabha scattering in the limit of fixed angle and high energy 
at the two-loop level. The integrals have been calculated using direct expansions 
of Mellin-Barnes representations, followed by a resummation of resulting harmonic 
series. 



1 Introduction 



The Bhabha scattering process, i.e. the elastic scattering of an electron and a positron is an 
invaluable tool in the luminosity determination in various experiments. Devices like the Interna- 
tional Linear Collider require a theoretical prediction for the cross section that contains, besides 
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exponentiated factors coming from the leading logarithmic terms, also fixed order contributions, 
and in particular the two-loop QED corrections. 

As far as large angle Bhabha scattering is concerned a prediction to leading order in the small 
electron mass (we shall later refer generally to any fermion species) has been derived in [1], 
although the logarithmic terms have been known before [2]. This prediction has been obtained 
from the massless result [3] by a careful matching procedure. Since heavy use has been made 
of results from different sources, it is definitely desirable to obtain the cross section by an 
independent method. Even though much progress has been made on the way to a result with 
an exact mass dependence [4,6,5], we are still far from the goal. In the mean time, we have 
found it much easier to evaluate the required diagrams in the same approximation as in [1]. 
The present paper presents our results for the set of planar Master Integrals (MI). The main 
tool are expansions of Mellin-Barnes (MB) integrals. 

In the next section we present some considerations on the choice of the Mis, we then describe 
the derivation and expansion of MB integrals. We subsequently give our results and finish with 
conclusions. An appendix contains our new MB representations for the five-line integrals. 



2 Master Integrals 



As explained in the introduction, it is our purpose to calculate the integrals in the limit of 
small fermion mass, which corresponds to high-energy scattering at fixed angle. Since this 
is a Sudakov approximation, we have to expect an overlap of collinear and on-shell infrared 
singularities leading to logarithms of the fermion mass reaching power four at the two-loop 
level. Even though the final result for the cross section will only be given at leading order, 
i.e. neglecting terms suppressed by m^, it is not immediately clear how many terms of the 
expansion of the individual Master Integrals are needed. Indeed, it is possible that some Mis 
appear with such coefficients that terms of relative order l/ni^ would result. These would, of 
course, have to cancel at the end but they would also make a higher order expansion in m^ 
necessary. One way to be sure that the choice of masters avoids this situation is simply to 
calculate the amplitude and find a set of Mis, whose coefficients would all be of relative order 
0(1) in the mass expansion. An alternative strategy, which can be applied without knowing the 
full amplitude is to require that the integrals have their mass dimension given by the kinematic 
invariants alone. To be specific, we shall require that for any integral J, the following relation 
holds: 

J = (m2)-2^s"F (in ( -— V -") + O (m^) , n = - dim / + 2e, (1) 



where s and t are the usual Mandelstam kinematic invariants, e is defined through the dimension 
of space-time in dimensional regularization, d = 4 — 2e, and F is some function that we wish 
to determine. The factor (m^)~^^ could have been replaced by (— s)^^^ through a reshuffling of 



logarithms, but we shall keep it as it is, because at intermediate stages of our calculation the 
mass has been set to unity for simplification. 

Whether Eq. (1) can be satisfied for a given integral or not depends entirely on the strength 
of the IR singularities. Indeed, usual Feynman graphs have only logarithmic branching points, 
but putting a dot on some line may result in stronger divergences. This is easily seen in a 
graph that would have just one massive line regulating the collinear divergence. Putting a 
dot on it, would be equivalent to a derivative, which transforms a logarithm ln"'(— m^/s) into 
nln"~^(— m^/s)/m^, which is precisely what we want to avoid. Of course, choosing Mis with 
property Eq. (1) is not yet a proof that there will be no coefficients with 1/m^, but we will be 
satisfied with the assumption, based on experience, that this is indeed the case. 

The consequence of the above considerations is that we have to slightly modify the set of 
integrals with respect to the original one presented in [5]. In particular, most of the dotted 
masters have to be replaced by those with irreducible numerators. The new, equivalent set, is 
given in Section 4. We have, of course, checked that the listed integrals are indeed algebraically 
independent. In case the reader would like to move back to the old set, we provide the transition 
formulae in [7]. 

Having chosen the set of integrals to calculate, there is another question to be considered, 
namely whether one needs to derive the expansions of the integrals both in the s and in the t 
channel, as they both occur in the Bhabha scattering amplitude. Fortunately, the limit under 
study is s,t ^ m'^, which means that the results are valid in both cases. For consistency of 
notation, one just has to change 

(2) 





and 

' ' 1/f-). (3) 



s t \s 

The above transformations imply some algebra in transforming the arguments of the polylog- 
arithms to a unique form. In order to spare the reader the effort, we provide the expansions in 
both channels in [7]. 



3 Expansion of Mellin-Barnes representations 



Mellin-Barnes representations for massive box integrals with the highest number of lines (seven) 
at the two-loop level have been derived in [8,9]. Since general powers of the propagators have 
been kept, it should, in principle, be possible to obtain representations for any box integral with 
a smaller number of lines. The occurrence of factors l/r(a), where a is a propagator power. 



seems to make some of the MB integrals trivially vanish. This is actually not the case, since 
they turn out to be singular in a and an analytic continuation in the relevant propagator power 
leads to a non-zero result with a smaller dimension of the integral. This procedure is not always 
optimal and it turns out to be better to derive MB integrals of smaller dimensionality directly 
for the considered graph. This is particularly true for the 5-line integrals. 

There is another case, where direct derivation of representations is necessary. We have seen in 
the previous section that we need master integrals with irreducible numerators. In such cases it 
is impossible to obtain the MB integrals directly from the results for the 7-liners. In any case, 
we derived our own representations by the standard technique of loop-wise integration, with 
Feynman parameter representations for tensors as given for example in [5]. The new integrals 
can be found in the Appendix. 

Having a MB integral at hand, one has to perform an analytic continuation in e from a range 
where the integral is regular to the vicinity of the origin, uncovering the singularity structure 
on the way. We perform this operation with the MATHEMATICA package MB [10]. As a result, we 
obtain multidimensional integrals, which still have a nontrivial dependence on both Mandelstam 
variables. At this stage, the integral can be cast into the following form 



(4) 




where the / function contains, amongst others, a product of F, or possibly ip functions, which 
have poles in z. Let us stress once more that, bar trivial cases, the / function is given by a 
multidimensional integral. Since it is difficult to directly take residues in this form, we change 
the order of integration and close the z contour to the right. This procedure is subsequently 
applied recursively, until no further poles at the required order of expansion occur. Since the MB 
package can perform numerical integrations, we could check some of the integrals beyond the 
leading order against the sector decomposition method with relatively large fermion masses. 
Moreover, we could confirm explicitly that the expansion is given by powers of rn^ and not of 
m, even though odd powers of the latter are found in individual terms of the series. 

Interestingly, the nontrivial dependence on the t/ s ratio occurred at the end only in one- 
dimensional integrals. This is in strong contrast to the massless calculation in [11], for example. 
The one-dimensional integrals were all of the harmonic type and could be done with XSUMMER 
[12]. The remaining, constant, integrals contained F functions of doubled argument, F(a±2z), 
which is a trace of the original mass. These have been dealt with high precision numerical 
integration followed by the use of the PSLQ algorithm [13]. The transcendental constants that 
could occur in the final result had to be Riemann C, numbers up to weight four, because of the 
correspondence between the expanded massive and the purely massless result. 



4 Results 



All our integrals are defined with the integration measure 
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In the presentation of the results, we will use the following notation 



which guarantees that our results are explicitly real in the Euclidean domain, s, t < 0; and 

SnA^) = ) \,. , / - In-^t) ln^(l - zt), Li„(^) = S„.i,i(^), (7) 

(n — ly.pl J t 



where Sn,p{^) i^ ^^^ Nielsen polylogarithm. We will, furthermore, always neglect the (tti^)"^'^ 
factor in front of the integrals, and keep the u Mandelstam variable in the prefactors, with the 
understanding that 

s + t + u = 0, (8) 



which is correct in the limit under consideration. In fact, we believe that the difficulties encoun- 
tered in the exact evaluation of the 5-line masters in our previous works are precisely connected 
to the occurrence of the u variable in the prefactors, which breaks the simplifications brought 
by the use of the conformally mapped variables 



^'l - Am^/s - 1 ^1 - 4mVt - 1 

'l - Am^/s + 1 ' Jl - Arn^/t + 1 ' 



^.1 7-line integrals 



The two 7-line topologies are given in Fig. 1. Here and in the following we will give the irreducible 
numerator in parentheses after the integral name. For the 7- and 6-liners we introduce the 
shorthand 

N={k2+p^f, (10) 
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Fig. 1. The planar 6- and 7-line topologies. The momentum distribution is common to all integrals. 

for the irreducible numerator as in [8,9]. For the calculation we used the MB representation 
from the mentioned papers. 



1 2L^ 
B714ml = + -—- 
e^ sH 
1 1 



1 



QsH 



"-10 L^ + 12 1} ln(x) + 6 L C2 + 6 Cs 

12 L^ + 28 L^ ln(x) - 4 L^ (-30 C2 + 3 \yi^{x)) 



- 4 L (9 Cs + 30 C2 ln(x) + \n\x) - 18 C2 ln(l + x) 

- 3 \n^{x) ln(l + a;) - 6 ln(a;) Li2(-a;) + 6 Li3(-a;)) 

-15C4-24C3ln(a;; 
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1 1 

e s^ 
1 



24s2 



13/6 L^ - 2 L^ ln(a;) - 2 L C2 - Cs 
29 L^ - 48 L^ ln(a;) - 408 L^ C2 



+ L [168 Cs + 384 C2 ln(x) + 16 \n^{x) - 288 C2 ln(l + x) 

- 48 \n^{x) ln(l + a;) - 96 ln(a;) Li2(-a;) + 96 Li3(-x)] + 144 (4 
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6sH 



-\} + L ln(a;) 

-2 L^ + 3 \} ln(a;) - L ln(x)2 - C2 ln(a;) 

-10 \} + 26 L^ ln(a;) - 2 L^ (-21 (2 + 9 ln^(a;)) 



- 2 L [30 Cs + 21 C2 ln(a;) - ln^(a;) - 18 C2 ln(l + x) 

- 3 \x\.{xf ln(l + a;) - 6 ln(a;) Li2(-a;) + 6 Li3(-a;)] 

+ 279 C4 + 18 Cs ln(a:) + 6 C2 ln'(a;) 
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B714m2(N) 
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e^ s 
1 1 
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e s^ 
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24^ 



L^ - L ln(a;) 



2 L^ - 3 L^ ln(a;) + L ln2(a;) + C2 ln(x) 



25 L^ - 48 L^ ln(a;) - 456 L^ C2 + L [192 Cs + 696 C2 ln(a;) 



+ 32 ln^(a;) - 288 C2 ln(l + a;) - 48 \\?{x) ln(l + x) 

- 96 ln(x) Li2(-a;) + 96 Li3(-x)] 

- 804 C4 - 168 Cs ln(a;) - 216 C2 \^{x) - 8 ln^(a;) 

+ 96 C2 ln(a;) ln(l + x) + 16 ln^(a;) ln(l + x) - 192 C2 Li2(-a;) 

+ 96 \w{x) Li3(-a;) - 192 Li4(-a;)| 
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e s 
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120i 

2 



-\} + L ln(x) 

-2 L^ a; + L^ [3 a; ln(x) - 1/2] - L a; ln2(x) - 4 C2 - x C2 ln(a;) 

- 125 L^ a; - 5 L"'^ [16 + 56 a; - 48 a; ln(a;)] 



5 \} [24 - 456 a; C2 - 144 x ln(a;)] 

5 L [120 C2 + 384 a; C2 + 192 a; Cs + 696 x (2 H^) 



96 X In (x) 



-32 a: ln(a;)2 - 288 x C2 ln(l 



a; 



48 X \n^{x) ln(l 



a; 



840 Cs 



- 96 a; ln(a;) Li2(-a;) + 96 a; Li3(-a;)] - 960 C2 - 1920 x C2 
+ 480 X Cs + 4020 a; C4 + 1440 x C2 ln(a;) + 840 x Cs ln(a;) 

+ 1080 a; C2 ln^(a;) + 80 a; ln-'^(a;) + 40 a; \n^{x) - 1440 a; C2 ln(l + 

- 480 X C2 ln(a;) ln(l + x) - 240 x In^(x) ln(l + x) 

- 80 X ln(x)^ ln(l + x) + 960 x C2 Li2(-x) - 480 x ln(x) Li2(-x) 

+ 480 X Li3(-x) - 480 x ln(x) Li3(-x) + 960 x Li4(-x)| 
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L^ - L ln(a;) 



2V x' + V (-1/4 + x/2 - 3 x^ ln(a;)) + L{l-x + x^ \n'{x)) 



-8 + 16 a; + 84 x^ - 48 x^ ln(x)) 



- 2 C2 + 4 a; C2 + a:^ C2 ln(x) 

+ ^|25L^x2 + L 
24 I 

+ L=^ (24 - 72 a;2 - 456 x^ (2 - 216 x"^ ln(x)) 

+ L (120 - 120 X - 60 C2 + 120 x (2 + 576 x^ (2 + 192 x^ (a 

+ 96 x^ \n{x) + 696 x^ C2 ln(a;) + 144 x^ \n^{x) + 32 x^ In^(x) 

- 288 x^ C2 ln(l + x) - 48 x^ ln(x)^ ln(l + x) 

- 96 x^ ln(x) Li2(-x) + 96 x^ Li3(-x)) 

+ |(3 ((-6 + 22 X + 60 x^) C2 + (-7 + 2 (7 - 6 x) x) Cs - 67 x^ (4) 

+ 2 x2 (-(ln(x) (54 C2 + 21 Cs + ln(x) (27 C2 + ln(x) (3 + ln(x))))) 
+ (9 + 2 ln(x)) (6 C2 + In'(x)) ln(l + x)))/6 
-2x2((4C2-3 ln(x)) Li2(-x) 

+ (3 - 2 ln(x)) Li3(-x) + 4 Li4(-x)) 
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4.2 6-line integrals 



The 6-line integrals are obtained from the 7-line topologies of Fig. 1, by removing line 1 for 
B613ml and line 3 for B613m2, and keeping the same momentum distribution. This follows from 
the fact that we use the same original MB representation from [8,9]. 
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1 1 

e 2si 
1 

24st 



L^ - 2L^ ln(x) + L[8C2 + ln^(a:)] 
7 L^ - 32 L^ ln(x) + L^ (24 C2 + 42 \n\x)) 



+ L (24 Cs - 72 C2 ln(x) - 16 \n%x) + 72 C2 ln(l + x) 
+ 12 In^(x) ln(l + x) + 24 ln(x) Li2(-x) - 24 Li3(-x)) 

/-672 C4 - 48 C2 ln^(a:) - In^(x) - 12 (6 C2 + ln2(x)) Li2(-x) 

72 Li4(-x) 



24stl 
+ 48 ln(x) Li3(-x 



;i7) 



B613ml(N) 



-L' - 4L + 2C2 



1 1 

e 2^ 

1/l3 + L2(-6-12 ln(a;)) 
OS I 



+ L (-24 + 24 C2 + 6 ln2(a;)) + 12 (2 + 12 Ca 



B613m2 = + ^{7 L^ - 28 L=^ ln(a;) + 6 L^ (26 (2 + 7 111^(3;)) 

+ L (-312 C2 ln(a;) - 28 ln^(a;) + 12 (6 C2 + ln^(a;)) ln(l + x) 

+ 24 \n{x) U2{-x) - 24 Li3(-a;)) 

+ 510 C4 + 156 C2 iri\x) + 7 \n\x) - 8 \n{x) (6 C2 + In^(x)) ln(l + x) 

+ 4 (6 C2 - 3 \n\x)) U2{-x) + 24 Li4(-x) 



(19) 



B613m2(N) 



1 1 

e 2s 
1 



-V 



K2 



^ 3 L-^ + 3 V (-2 - 6 ln(a;)) + 3 L (18 C2 + 6 \n\x)) 

6s L 

- -|-9 Cs - 10 C2 ln(x) - \n\x) + 6 C2 ln(l + x) + In^(x) ln(l + x) 
+ 2 ln(a;) U2{-x) - 2 Li3(-a;) 



(20) 



4.3 5-line integrals 



The topologies for the 5-hners with numerators are given in Fig. 2. For the evaluation of the 
integrals we used our own MB representations given in the Appendix. The position of the dots 
on the dotted integrals (having a "d" in their name) has been defined in [5]. These integrals 
belong to our original MI set. 



1 1 rri 

B512ml = + — - [L] 



e^ s 
1 1 

e s 
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6s 



-3/2 L^ + L (-2 + ln(x)) + C2 
-7 L^ + V (-18 + 9 ln(a;)) + L (-24 + 24 C2 + 12 ln(x) - 3 ln^(a;)) 



+ 12 C2 + 12 Cs - 24 C2 ln(x) - ln3(a;) + 18 C2 ln(l 



X) 



+ 3 In^ix) ln(l + x) + 6 \n{x) Li2(-a;) - 6 Li3( 
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Fig. 2. The 5-line topologies. The momentum distribution has been chosen to make the derivation of 
MB representations easier. 
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-3 L^ + 6 L ln(x) - 3 In^(x) - 24 C2 

(-4 L^ - L^ (6 - 9 ln(x)) - L (6 C2 - 12 ln(x) + 6 \n^{x)) 

- 48 C2 - 54 Cs - 18 C2 ln(x) - 6 In^(x) + In^(x) 
+ 18 C2 ln(l + x) + 3 In^(x) ln(l + x) 

+ 6 ln(a;)Li2(-x)-6Li3(-s)| 
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B512m2(k2 ■ Ps) = H 



3L2 + L (12 + 6X-6 ln(x)) 



+ 24 C2 - 12 ln(x) - 6 a; ln(a;) + 3 ln2(x) 

+ -^{4 L^ + L^ (3 (8 + 3 x) - 9 ln(x)) 

+ L [6 (10 + 5 x + C2) - 12 (3 + x) ln(x) + 6 In^(x)] 

+ 12 C2 + 54 C3 - 60 ln(x) - 30 x ln(x) + 18 C2 ln(x) + 12 \n^{x) 

+ 6 X \n^{x) - In^(x) - 18 C2 ln(l + x) - 3 \n^{x) ln(l + x) 

-6 ln(x)Li2(-x) + 6Li3(-x)| 



(23) 



B512m2d2 



: + — (-12 L^ + 9 L^ ln(x) - 36 L C2 
Qst L 

- 12 Cs - 12 C2 ln(x) - ln=^(x) + 18 C2 ln(l + x) 

+ 3 In^(x) ln(l + x) + 6 ln(x) Li2(-x) - 6 Li3(-x)| 
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B512m3 = + 



12m I 



-6L2(6C2 + ln'(x)) 



10 



- 6 L (-4 Cs + 4 C2 ln(a;) - 12 (2 ln(l + x) - 2 In^ix) ln(l + x) 

- 4 \n{x) Li2(-x) + 4 Li3(-a;)) 

+ 312 C4 + 72 Cs Hx) + 36 C2 \n\x) + \n\x) - 24 (3 ln(l + x) 
+ 24 C2 ln(a;) ln(l + x) - 36 C2 ln^(l + a;) - 6 ln^(a;) ln^(l + x) 



24 ln(x) Si,2(-x) + 12 (8 C2 + In'(x) 
48 ln(a;) Li3(-a;) + 24 ln(l + x) Li3( 

72 Li4(-a:) + 24 82^ 



2 ln(a;) ln(l + x)) Lial 
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-10 L^ 



2 ' 2 L ln(a;) - C2 
12 L^ \n{x) + 
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- 21 Cs - 42 C2 ln(a;) - 2 In^(x) + 36 C2 ln(l + x) 

+ 6 ln=^(a;) ln(l + x) + 12 ln(x) Li2(-x) - 12 Li3(-x)| 

+ e — (-5 L^ + 10 L^ ln(x) + 111 L^ (2 
3st I 

+ L (104 Cs - 126 C2 ln(x) - 6 In^(x) + 108 C2 ln(l + x) 
+ 18 In^(x) ln(l + x) + 36 ln(x) Li2(-x) - 36 Li3(-x)) 

- 372 C4 - 78 Cs ln(x) + 30 C2 ln^(a;) + 2 In^(x) + 12 Cs ln(l + x) 

- 60 C2 ln(x) ln(l + x) - 8 In^(x) ln(l + x) + 18 C2 ln^(l + x) 
+ 3 ln2(x) ln=^(l + x) + 12 ln(x) Si,2(-x) 

- 6 (4 C2 + 3 \n\x) - 2 ln(x) ln(l + x)) Li2(-x) 
+ 24 ln(x) Li3(-x) - 12 ln(l + x) Li3(-x) 

-12Li4(-x) -12S2 2(-a:' 
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B512m3(k2 • P2) = t ( - ) j L^ (6 x C2 + 2 x ln(x) + 2 x^ ln(x) + x In^(x)) 
+ L (16 X C2 - 8 x^ C2 - 4 X Cs - 2 ln(x) + 2 x^ ln(x) 



4 X C2 ln(x) + 2 X In^(x) - 2 x^ In^(x) - 12 x C2 ln(l 
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-2x In^(x) ln(l 
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4 X ln(x) Li2( 
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4x Li3 
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+ T7^ - +120 C2 + 360 X C2 - 120 x' C2 - 1560 x C4 - 480 x Cs 
120 \u/ I 

- 240 x^ Cs - 240 X C2 ln(x) - 480 x^ C2 ln(x) - 360 x Cs ln(x) 
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+ 30 In^ix) + 60 X \b.^{x) - 30 x^ \n^{x) - 180 x (2 ln^(a;) 

- 20 a; ln^(a;) - 20 x^ \n\x) -5x \n^{x) + 720 x^ C2 ln(l + x) 

+ 120 X Cs ln(l + a;) - 120 x C2 ln(a;) ln(l + a;) + 120 x^ ln^(a;) ln(l + x) 
+ 180 a; C2 ln^(l + a;) + 30 x ln2(x) 111^(1 + x) + 120 x ln(x) Si,2(-x) 
+ 60 X (-8 C2 - ln^(a:) + 2 ln(x) (2 x + ln(l + x))) Li2(-x) 

- 240 x^ Li3(-x) + 240 x ln(x) Li3(-x) - 120 x ln(l + x) Li3(-x) 

-360xLi4(-x) -120xS2,2(-a;)j (27) 
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+ x^ ln2(x)) + L (-4 C2 + 20 X C2 - 8 X Cs - 4 x^ Cs 

- 2 ln(x) + 2 x^ ln(x) + 8 x C2 ln(x) + 4 x^ C2 ln(x) 
+ 4 X In^(x) - 24 X C2 ln(l + x) - 12 x^ C2 ln(l + x) 

- 4 X In^(x) ln(l + x) - 2 x^ In^(x) ln(l + x) 

- 4 X (2 + x) ln(x) Li2(-x) + 4 x (2 + x) Li3(-x)) 

+ ^ (-] |-24 (-1 + (-3 + x) x) C2 + 72 Cs + 24 X Cs - 552 x (2 + x) C4 

- 48 X (2 + x) Cs ln(x) + ln(x) (48 (1 + x) C2 

+ ln(x) (6 - 6 X (-2 + X + 16 C2 + 8 X C2) + 4 (1 + x) ln(x) 

- X (2 + x) ln2(x))) + 24 X (2 + x) Cs ln(l + x) 

+ 4 (-3 (1 + 3 x) (6 C2 + In'(x)) + In(-x) (24 x (2 + x) C2 

+ 3 ln(x) (2 + 6 X + X (2 + x) ln(x)))) ln(l + x) 

+ 6 X (2 + x) (6 C2 - 2 In(-x) ln(x) + ln2(x)) ln^(l + x) 

+ 96 X (2 + x) C2 (C2 - In(-x) ln(l + x) - Li2(-x)) 

+ 24 ln(x) (C2 - In(-x) ln(l + x) - Li2(-x)) 

+ 72 X ln(x) (C2 - In(-x) ln(l + x) - Li2(-x)) 

+ 12 X (2 + x) \n^{x) (C2 - In(-x) ln(l + x) - Li2(-x)) 

+ 24 Li3(-x) + 72 X Li3(-x) + 48 x (2 + x) ln(x) Li3(-x) 

- 24 X (2 + x) ln(l + x) Li3(-x) - 12 x (2 + x) (2 Cs ln(x) 

+ 2 C2 ln(x) ln(l + x) - In(-x) ln(x) ln^(l + x) - 2 ln(x) Si,2(-x) 

- 2 ln(x) ln(l + x) Li2(-x) + 6 Lu{-x) + 2 S2,2(-a;))| (28) 



B513m = + — (l^ (-18 C2 - 3 ln2(x)) + L (12 Cs + 2 In^(x) + 36 C2 ln(l + x) 
+ 6 lii^ix) ln(l + x) + 12 ln(x) Li2(-x) - 12 Li3(-x)) 
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B513md3 



- 81 C4 - 12 Cs ln(l +x) -2 In^(x) ln(l + x) - 18 C2 ln^(l + x) 

-3 ln^(a;) 111^(1 +a;) - 12 ln(a;) Si,2(-a;) -6 ln(a;) (ln(a;) 

+ 2 ln(l + x)) U2{-x) + 12 ln(x) Li3(-a;) + 12 ln(l + x) Li3(-a;) 

-12Li4(-a:) + 12S2,2(-a;)| 
1 1 



+ 



e st 
3x I 



-V + L ln(a;) 



L^ + 12 L^ \n{x) + L (24 C2 - 3 \n'^{x)) 



- 21 C2 ln(x) - In^(x) + 18 C2 ln(l + x) + 3 \n^{x) ln(l + x) 
+ 6 ln(a;) Li2(— a;) — 6 Li3(— x) 



(29) 



(30) 



B513m(ki • P2) = + — - 



1 /s 
12 \u. 

+ L^ (-18 a; C2 + 6 \n{x) + 6x \n{x) -3x Inixf) 
+ L (36 C2 - 36 X C2 + 12 X Cs + 6 ln(a;) - Q x"^ \n{x) 

- 12 a; ln(a;)^ + 2 x ln(a;)^ + 36 a; C2 ln(l + a;) + 6 a; ln(x)^ ln(l + x) 
+ 12 a; \n{x) Li2(-a;) - 12 x Li3(-a;)) + 12 (2 - 12 x C2 

+ 12 x^ C2 - 12 Cs + 12 a; Cs - 81 X C4 - 3 ln(x)2 - 6 x ln(x)^ 

+ 3 x^ ln(x)^ - 2 ln(x)^ + 2 x ln(x)^ - 36 C2 ln(l + x) 

+ 36 X C2 ln(l + x) - 12 X Cs ln(l + x) - 6 \n{xf ln(l + x) 

+ 6 X ln(x)^ ln(l + x) - 2 x ln(x)^ ln(l + x) - 18 x C2 ln(l + x)^ 

— 3 X ln(x)^ ln(l + x)^ — 12 x ln(x) Si^2(— 2;) — 6 ln(x) (2 — 2 x + x ln(x) 
+ 2 X ln(l + x)) Li2(-x) + 12 Li3(-x) - 12 x Li3(-x) + 12 x ln(x) Li3(-x) 

+ 12 X ln(l + x) Li3(-x) - 12 x Li4(-x) + 12 x S2,2(-a;) j (31) 



B513m(ki • Pa) 



^ 12 U 
+ L^ (18 x^ C2 - 6 X ln(x) - 6 x^ ln(x) + 3 x^ ln(x)2) 

+ L (-72 X C2 - 12 x^ Cs + 6 ln(x) - 6 x^ ln(x) - 6 x ln(x)^ 

+ 6 x^ ln(x)^ - 2 x^ ln(x)^ - 36 x^ C2 ln(l + x) - 6 x^ ln(x)^ ln(l + x) 

- 12 x^ ln(x) Li2(-x) + 12 x^ Li3(-x)) - 12 C2 - 60 x C2 - 12 x^ C2 
+ 24 X Cs + 81 x^ C4 - 12 X C2 ln(x) - 12 x^ C2 ln(x) - 3 ln(x)^ 

- 6 X ln(x)^ + 3 x^ ln(x)^ + 2 x ln(x)^ - 2 x^ ln(x)^ 

+ 72 X C2 ln(l + x) + 12 x^ Cs ln(l + x) + 12 x ln(x)^ ln(l + x) 

+ 2 x^ ln(x)^ ln(l + x) + 18 x^ C2 ln(l + x)^ + 3 x^ ln(x)2 ln(l + xf 
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+ 12 x^ \n{x) Si,2(-a;) + 6 x ln(x) (4 + x ln(a;) +2x ln(l + x)) Lisi 



-X) 



- 24 X Li3(-x) - 12 x^ ln(x) Li3( 



-X 



12x2 ln(l + x) Li3(-x) 



+ 12 x^ Li4( 



-X 



12x^8 



2,21 



-X 



(32) 



B513m(k2 • ps) = + — - 



-(- 

12 Vm 



+ L^ (18 X C2 + 6 X ln(x) + 6 x^ ln(x) + 3 x \n^{x)) 

+ L (36 X C2 - 36 x^ C2 - 12 X Cs - 6 ln(x) + 6 x^ ln(x) - 12 x^ hi^{x) 

— 2 X In (x) — 36 X C2 ln(l + x) — 6 x In (x) ln(l + x) 

- 12 X ln(x) Li2(-x) + 12 x Li3(-x)) - 12 C2 + 12 x (2 

- 12 x^ (2 - 12 X Cs + 12 x^ C3 + 81 X C4 + 3 ln=^(x) 

+ 6 X ln2(x) - 3 x^ ln2(x) - 2 x In^(x) + 2 x^ In^(x) - 36 x (2 ln(l + x) 
+ 36 x^ C2 ln(l + x) + 12 X Cs ln(l + x) - 6 x \n^{x) ln(l + x) 
+ 6 x^ \n^{x) ln(l + x) + 2 x In^(x) ln(l + x) + 18 x C2 ^y?{1 + x) 
+ 3 X In (x) In (1 + x) + 12 x ln(x) Si,2(— a;) + 6 x ln(x) (ln(x) 
+ 2 (-1 + X + ln(l + x))) Li2(-x) + 12 x Li3(-x) - 12 x^ Li3(-x) 

— 12 X ln(x) Li3(— x) — 12 x ln(l + x) Li3(— x) 

+ 12xLi4(-x) -12xS2,2(-a;)j (33) 



B514m 



1 1 

e^ s 
1 1 

e s 



[L] 



-3/2 L^ + L (-2 + ln(x)) + (2 



- — (-7 L^ + 1} (-18 + 9 ln(x)) + L (-24 + 24 (2 + 12 ln(x) - 3 Xn^ix)) 
6s I 

+ 12 C2 + 12 Cs - 24 C2 ln(x) - \n\x) + 18 C2 ln(l + x) 
+ 3 ln2(x) ln(l + x) + 6 ln(x) Li2(-x) - 6 Li3(-x)| 



(34) 



B514md = 1-11 1^ + 211^ ln(x) + L (30 C2 - 9 \n^{x)) 

- 24 C2 ln(x) - ln2(x) + 18 C2 ln(l + x) 
+ 6 ln(x) Li2(— x) — 6 Li3(— x) 



3 \n^{x) ln(l 



X 



(35) 
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5 Conclusions 



With tlie lielp of a direct expansion of MB integrals, we obtained all of the planar box masters to 
leading order in the expansion in the small fermion mass. The calculation has been substantially 
simplified by the fact that all nontrivial integrations were one-dimensional. Unfortunately, the 
remaining, non-planar Mis introduce two further complications. First, the final integrals are 
not directly one-dimensional. This is per se not an insurmountable obstacle, since a similar 
situation has been dealt with during the evaluation of the fully massless masters. However, in 
our case it seems non-trivial to bring the integrals to the form of harmonic sums. The second 
problem is connected to the fact, that MB representations for the non-planar case involve all 
three kinematic invariants. All in all, despite the above problems we hope that the technique 
used in this paper will prove successful also in this case, and we will be able to provide the 
complete Bhabha scattering cross section in the limit of high energy and fixed angle from 
a direct diagrammatic calculation and thus verify results obtained by matching between the 
massive and massless cases [1]. Moreover, our approach can be applied to problems with more 
scales, in particular, we have applied it to the integrals occuring in the Bhabha cross section 
with two fermion species [14]. 
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A Mellin-Barnes representations for 5-line integrals 



This appendix contains the MB representations that we directly derived for the 5-line integrals 
with as well without irreducible numerators. Notice that the momentum p^ can be any of the 
external momenta, because we performed the integration with a single loop momentum in the 
numerator with vector index kept free. As far as the notation is concerned, the Oj denote the 
power of line i. The numbering can be read off Fig. 2. 



B512m2 = ^ — FT^TFT^ 7; TTTTZ^ I da j d(3 j d-f (A.l^ 



UU^[a^]T[4-2e-a^3]{27^^)^ 
r[-7]r[2 - e - a245 - 7 - tt]r[-a] 
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T[a2 + a]T[a4 + a]r[4 — 2e — ana — P] 

r[-2 + e + 0245 + 7 + a - /3]r[ai + /5] 

r[4 - 26 - 02245 -2a + P]T[2 - e - 024 - 7 - « + /3] 

r[4 - 2e - 0245 + /^]r[4 - 2e - 022445 - 2-f - 2a + p] 



B512m2(pe • ka) 



' _ ]^ ')'*12345 -Ig2e7£; 



r[oi]r[o2]r[o3]r[o4]r[o5]r[4 - 2e - oi3](2 



TTi 



+ioo +ioo +ioo 

da I dp d'j 



— s 



l(2-e-a245-a+/3-7) 



-trr[-a]r[-7]r[oi + /3] 



r[4 - 2e - 0113 - P]T[-2 + e + ai3 + P] 
r[5 - 2e - 0245 + /9] 
r[-2 + e + 0245 + a - /9 + 7] r[2 - e - O245 - « - 7] 

r[4 - 2e - 022445 - 2a + /5 - 27] r[5 - 2e - 022445 - 2a + /3 - 27] 

Pe ■ {Pi + P2)r[5 - 2e - 022445 + /3 - 27] 

r[4-2e-022445-2a + /5-27] 

r[o2 + 7]r[o4 + 7]r[3 - e - 024 - « + /? - 7] 

+ r[4 - 2e - 0224445 + /3 - 27]r[5 - 2e - 022445 - 2a + /3 - 27] 

— rtn^ A- -I- /9 ^1 



r[2 - e - 024 - « + /3 - 7] 

• Pi r[i + 02 + 7] r[o4 + 7] + Pe • P3 r[o2 + 7] r[i + 04 + 7] 



p. 



(A.2) 



B512m3(pe ■ ks) 



'_]^')«12345 g2e7B 



n-=irNr[4-2e-Oi23](27r2)4 



lOO +100 

(ia / dp 



(^7 / (i5 



-^) 



(4-2e-ai2345 -/3-5-7) 



r[-/5] r[-7] r[-5] r[o3 + a + 2 /?] r[2 - e - 045 + « - 5 - 7] 



r[2 



ai3 



r[7-3e-oi2345 - P] 
/?]r[2-e-023-a-/9] 



r[o5 - a + 2 7] r[l + 05 - a + 2 7] 
r[-4 + 2e + oi2345+/3 + 5 + 7] 

|r[4-2e-oi235-/9-5-7] 

(pe ■ P2) r[i + 05 + 7] r[-a + 7] - (pe • Pi) r[o5 + 7] r[i - a + 7] 

r[o5 - a + 2 7] r[l + 05 - a + 2 5 + 2 7] 

+ [(Pe • Ps) - (Pe ■ Pi)] r[5 - 2e - 01235 - /3 - (5 - 7] 

r[o5 + 7] r[-a + 7] r[l + 05 - a + 2 7] r[o5 - a + 2 (5 + 7)]! (A.3) 
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. r, +«oo +«oo +ioo +ioo 

B513m(pe • ks) = —^ — ^^ \^^' — TTTT" ^ j da j d(3 j d-f j d6 



J ' ' ' ' _ . . „jQQ — lOO — JOO — lOO 



n,tirHr[4-2e-ai23](27r2)4^ 

r[2 - e - ai3 - a - 7]r[ai + 7] 

r[-4 + 2e + ai2345 + a + /J + (5] 

r[6 - 3e - ai2345 - «] 

r[-«]r[-/5] 



r[4 - 2e - ai2345 -a-^-l3-5] 



r[7 - 3e - ai2345 - a] 

r[-^] 

r[4 - 2e - an23 - 2a - 7] 

r[2 — e — ai2 — a] 
r[8 - 4e - 0112233445 - 2a - 7 - 2/5 - 25] 

r[4 - 2e - aii23 - 7] 

r[9 - 4e - aii2233445 - 2a - 7 - 2/? - 25] 

{-{{Pe ■ Pi) r[7 - 3e - ai2345 - tt]r[8 - 4e - 0112233445 - 2a - 7 - 25] 

r[9 - 4e - aii2233445 - 2a - 7 - 2/5 - 25]r[4 - 2e - 01234 - a - /3 - 5] 

r[a4 + 5]r[-2 + e + 0123 + a + 7 + 5]) 

+r[6 - 3e - ai2345 - «] 

{{Pe ■ {Pi + P2)) r[9 - 4e - aii2233445 - 2a - 7 - 25] 

r[8 - 4e - aii2233445 - 2a - 7 - 2/5 - 25] 

r[5 — 2e — ai234 — a — /3 — 5]r[a4 + 5] 

r[-2 + e + ai23 + « + 7 + '^] 

+r[8 - 4e - aii2233445 - 2a - 7 - 25] 

r[9 - 4e - aii2233445 - 2a - 7 - 2/5 - 2 * 5] 

r[4 - 2e - ai234 - a - /5 - 5] 

((pe ■ P3) r[l + 04 + 5]r[-2 + e + ai23 + a + 7 + 5] 

+ (Pe ■ Pi) r[a4 + 5]r[-l + e + ai23 + a + 7 + 5]))) (A.4) 



, , D +Joo +Joo +Joo 



''''"(P--"^^ = n^irMr[5-2e-ai23](2vr.)^ J da J d(3 J d^ J d6 

/'_g\(4-2 e)-ai2345-a-/3-5 ('_^)'5 

r[-4 + 2 e + ai2345 + a + /5 + 5] 
r[6 -3 e-ai2345 - «] 

r[-a] r[-/5] 

r[7 - 3 e - ai2345 - a] r[5 - 2 e - 0123] 

n-5^ 

r[4 - 2 e - aii23 - 2 a - 7] r[5 - 2 e - 01123 - 2 a - 7] 
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r[2 - e - ai3 



a 



7] 



r[8 - 4 e - aii2233445 -2a-2p-26--f] 

r[4 - 2 e - ai2345 -a-P-S--f] 
r[9 - 4 e - aii2233445 -2a-2(3-26-^] 

I {Pe ■ Ps) r[l + 04 + 5] r[6 - 3 e - ai2345 - "] 

r[4 — 2 e — ai234 — a — /? — 5] r[3 — e — ai2 — a] 
r[8 - 4 e - aii2233445 -2a-26--f] 

r[9 - 4 e - 0112233445 - 2 « - 2 /3 - 2 5 - 7] 

r[5 - 2 e - aii23 - 7] r[4 - 2 e - 01123 - 2 a - 7] 
r[ai + 7] r[-2 + e + ai23 + a + 5 + 7] + T[a4 + S] 

-{pe ■ Pi) r[7 - 3 e - ai2345 - a] 

r[4 - 2 e - ai234 - « 

r[8 — 4 e — ail2233445 
r[9 — 4 e — ail2233445 

r[3 - e - ai2 - a] r[5 - 2 e - 01123 - 7] 

r[4-2e-aii23-2a-7]r[ai + 7] 
+r[2 - e - ai2 - a] r[4 - 2 e - 01123 - 7] 

r[5 - 2 e - aii23 - 2 a - 7] r[l + ai + 7] 

r[-2 + e + ai23 + a + 5 + 7] 

+r[6 - 3 e - ai2345 - a] r[3 - e - ai2 - a] 

r[5 - 2 e - aii23 - 7] r[4 - 2 e - 01123 - 2 a - 7] 



2 a - 2 5 - 7] 
2a-2p-26--f] 



r[ai+7] 



{{Pe ■ {pi + P2)) r[5 - 2 e - ai234 



a 



P-5] 



r[9 - 4 e - aii2233445 -2a-26--f] 

r[8 - 4 e - aii2233445 -2a-2p-26-^] 

r[-2 + e + ai23 + a + 5 + 7] 

+ {pe ■ pi) r[4 - 2 e - ai234 - a - (3 ~ 6] 

r[8 - 4 e - aii2233445 -2a-26--f] 

r[9 - 4 e - aii2233445 -2a-2P-25-^] 

r[-l + e + ai23 + a + (5 + 7] 



(A.5) 
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